Interfaces in driven Ising models: shear enhances confinement 
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We use a phase-separated driven two-dimensional Ising lattice gas to study fluid interfaces exposed 
to shear flow parallel to the interface. The interface is stabilized by two parallel walls with opposing 
surface fields and a driving field parallel to the walls is applied which (i) either acts locally at 
the walls or (ii) varies linearly with distance across the strip. Using computer simulations with 
Kawasaki dynamics, we find that the system reaches a steady state in which the magnetization 
profile is the same as that in equilibrium, but with a rescaled length implying a reduction of the 
interfacial width. An analogous effect was recently observed in sheared phase-separated colloidal 
dispersions. Pair correlation functions along the interface decay more rapidly with distance under 
drive than in equilibrium and for cases of weak drive can be rescaled to the equilibrium result. 
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Shearing a fluid is a fundamental example in the statis- 
tical mechanics of systems driven away from equilibrium 
PQ. Profound effects on the liquid structure [2 and new 
phase transitions [3 can result. In a recent intriguing ad- 
vance in colloidal science interfacial fluctuations of capil- 
lary wave- type were observed by direct optical means [4] ; 
these fluctuations can be strongly supressed by shearing 
in the direction parallel to the interface [5] . 

Kinetic Ising models which have been driven out of 
equilibrium are of continuing importance in condensed 
matter physics [6j [8] . Underlying is the idea of coarse 
graining in time [9] , where a record of the dynamical state 
is only kept at discrete points in time, which is an ap- 
proximation that can often be justified by the quality and 
range of deductions made. Mass is conserved, but there 
is no inertia and associated hydrodynamics. Neverthe- 
less currents are predicted as well as densities. Examples 
showing the value of such stochastic models include stud- 
ies of interfaces under spatially homogeneous and random 
driving fields [7], domain growth [8] and nucleation [10] 
under shear, as well as of the precursor film in complete 
wetting [TT] . 

In this Letter, we investigate the effects of shear on 
the interface between coexisting phases using the two- 
dimensional (2d) Ising lattice gas with Kawasaki dynam- 
ics [12]. The interface is induced by fixing spins on the 
boundaries of a strip; the analogue of shearing in a lattice 
gas will be explained below. Studying the 2d problem is 
interesting in its own right, because in equilibrium the in- 
terface is known to be rough for all temperatures T below 
the bulk critical temperature T c . If the interface is con- 
fined to a strip of width L yi then it sweeps out essentially 
the entire strip, meandering back and forth between the 
walls, such that the interfacial width w ~ L yi despite 



the entropic repulsion from the walls [13]. Consequently, 
at fixed T, the magnetization profile as a function of 
distance y across the strip scales as m\ ) (T)M. e(l (y / L y ) , 
where rat>(T) is the spontaneous magnetization in bulk 
and A^eq is the finite size scaling function for a strip with 
size L x (^> L y ) along the walls; "eq" labels equilibrium 
quantities. 

Using Monte Carlo (MC) simulations we find that the 
2d Ising model under shear-like drive parallel to the in- 
terface achieves a steady state that is characterized by 
an effective length scale L*, such that the magnetization 
profile m(y) at given T obeys 

^ Me « (l0 +M ™"^ with L *y < L y> W 

where M cor r(y) is a boundary correction term, which de- 
cays away from the boundaries on the scale of the bulk 
correlation length £ (on which the exponential decay of 
spin-spin correlations occurs in the homogeneous sys- 
tem). A conjecture based on Eq. is that shear-like 
drive acts as effective confinement of the system. 

To be specific, in spin language the Hamiltonian is 
H = —J^Z <i j > o'iO'j^ where the sum runs over nearest 
neighbour sites z, j; J > is the spin- spin coupling con- 
stant, and the spins take on values Oi = ±1, correspond- 
ing to particle occupation numbers ti = (cr^ + l)/2 = 0, 1. 
The interface is induced and localized by two walls of 
fixed spins <r^ = +1 at the top (y = L y + 1) and <r^ = — 1 
at the bottom (y = 0) edges of the strip. Periodic bound- 
ary conditions are applied in the ^-direction. In order to 
induce shear, we apply a force field JF(y) parallel to the 
x-direction. In model I only the particles in the layer 
adjacent to each wall are driven along the walls and in 
opposite directions: F(L y ) = Fq and F(l) = — Fq, and 
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F{y) — otherwise; for F — > oo these layers form an 
asymmetric exclusion processes 0, coupled to an Ising 
strip. In model II the force varies linearly across the strip, 
such that it vanishes in the middle of the slit: F(y) = 
w[y — (Ly + l)/2], where to = dF(y)/dy is the (dimension- 
less) field difference between adjacent rows; we also use a 
scaled variable y = (2y—L y — l)/L y , see Fig.JT^L. Model II 
can be thought of as mimicking the effects caused by the 
flow of a background solvent [5] . The system evolves un- 
der spin-exchange Kawasaki dynamics [12 j, correspond- 
ing to particle hopping to an (empty) nearest neighbour 
site. The work done by (or against) the force field, AF, 
in a trial move enters a modified Metropolis acceptance 
rate, min{l, exp(-(A# + AF)/(k B T))}, where AH is 
the change in internal energy and k B is the Boltzmann 
constant. The dynamics capture the local conservation of 
particle number and the competition of forced transport 
with diffusive motion. 

We have performed extensive simulations using single- 
spin and multi-spin coding techniques [14] , the latter gen- 
eralized to include drive; this facilitates simultaneously 
running 64 systems per processor core. We present here 
results for L x = 200 and L y varying from 10 to 40 at 
fixed total magnetization m = <Ji = for T/T c = 
0.75,0.85, and 0.95 (where T c = 2.2619 J/k B ), so that 
Ly/t; > 1. Sampling large-scale interfacial fluctuations 
with inherently slow Kawasaki dynamics is challenging 
and requires run lengths of the order of Nmc = 10 8 
MC steps (L x x L y trial moves form one MC step). For 
each run we perform initially TVmc MC steps and find 
that thereafter all observables of interest fluctuate around 
their mean value. We conclude that a steady state is 
reached, in which we gather statistics for a further Nmc 
steps. Statistical errors are obtained with blocking and 
bootstrap techniques. Fig. [I] displays snapshots of con- 
figurations [T5] . 

Driving creates a non- vanishing current profile j(y) 
parallel to the walls; j(y) gives the net (mean) num- 
ber of particles that move from x to x + 1 at given y 
in one MC step. For the boundary-driven case (model 
I) the current is localized at both walls, y = l,L y , and 
vanishes in the middle of the system (apart from a very 
small anti-current at y = 2, L y — 1). In the bulk-sheared 
case (model II) j(y) varies smoothly with y, see Fig. [l]i. 
For small values of oo a near-linear behaviour is observed. 
For strong drive the current saturates upon approaching 
the walls. Surprisingly \j(y)\ displays pronounced max- 
ima offset from the center of the strip, indicating an op- 
timum between sufficiently strong drive and sufficiently 
large number of particle- hole pairs. 

Having established the occurence of shear flow, we in- 
vestigate the interfacial magnetization profile m(y) = 
L~ 1 (Y JX cr(x,y)), where the angles denote an average in 
the steady state. We have checked that the aspect ra- 
tio L x /L y is large enough such that at equilibrium the 
exact result [13] for m(y), for L x — > oo, is indistinguish- 
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FIG. 1: (color online) a) Illustration of the 2d Ising strip be- 
tween +/— walls under an external driving field JF(y) paral- 
lel to the walls. The field acts with strength Fo at the bound- 
ary layers only (model I) or varies linearly with (scaled) dis- 
tance y across the strip, where oo determines the slope (model 
II). b) Snapshot from simulation for T/T c = 0.75, L x = 200, 
and L y = 20 at equilibrium; black (white) regions indicate 
(Ti — +1( — 1). c) Same as b) but under strong drive according 
to model II with uj = 2. d) The current profile j(y) parallel 
to the walls as a function of y for different values of oo (as in- 
dicated) in model II for the same values of T, L x , L y as in b. 



able within statistical errors from our simulation results 
(not displayed) using Glauber dynamics [16], where m 
can fluctuate. Using Kawasaki dynamics (that constrain 
m = 0) leads to larger flat regions near the walls and 
a sharper interfacial region, see Fig. |2^i; this effect in- 
creases strongly if L y /^ ^> 1. The importance of the 
contribution of capillary- wave type fluctuations to m(y) 
can be gauged by the stark contrast to the kink-like be- 
haviour obtained in equilibrium mean-field theory. For 
cases of weak boundary drive (small values of Fo) the 
magnetization near the walls increases and the interfacial 
region gets squeezed (becomes narrower) upon increasing 
Fo; this behaviour is more apparent for higher values of 
T. For Fq > 5 there are no further changes in the pro- 
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FIG. 2: (color online) Scaled magnetization profiles 
m(y)/mb(T) as a function of the scaled distance y for T /T c — 
0.75, L x = 200, L y = 20. a) Model I: Kawasaki equilib- 
rium result, Glauber equilibrium result, mean-field result, 
Kawasaki result for Fo = 50, and the latter rescaled as 
ra(0.83y). Inset: Variation of the interfacial width w with Fo 
for three different temperatures (as indicated), b) Kawasaki 
simulation results for model II are shown in equilibrium and 
for a range of values of the field gradient co (as indicated), as 
well as using rescaling m(a±y), where a± = 1/1.15, 1/3.4 for 
uo = 0.025, 1, respectively. Error bars are of the order of the 
line thickness. Inset: Variation of the scaled interfacial width 
w/Ly with the scaling variable l/(L y uo s ), where s = 0.3, for 
uo = 0.025 (diamonds), 0.1 (filled circles), 0.2 (squares), 0.25 
(stars), 0.5 (crosses), and 1 (triangles). 



file within error bars; in Fig. [2^i the limiting behaviour, 
Fq = 50, is shown. The effect on m(y) arises from a com- 
petition of advection and diffusion: the boundary drive 
helps to break up clusters into smaller constituents; these 
possess high mobility and are hence able to migrate to- 
wards and coalesce with their bulk phase for energetic 
reasons. This dynamic mechanism makes the observed 
additional confinement of the interface physically rea- 
sonable. The effect of the boundary drive cannot be 
mimicked by assuming equilibration in static boundary 



fields acting at y = 2 and L y — 1. Even infinitely strong 
surface fields h(y = 1) = —h(y = L y ) = oo induce a 
significantly weaker effect. Particularly striking is the 
scaling behavior of m(y), see Eq. ([!]), which indicates 
that the effect of the boundary drive on the interfacial 
region of the profile is the same as that induced by in- 
creasing the confinement, and hence reducing the width 
of the strip from L y to L* = 0.83L y for Fo > 5. This 
corresponds to the scaling m(a±y) « m eq (i/), where for 
strong boundary drive (such that essentially Fq — > oo) 
we find a± = Ly/L y = 0.83 for L x = 200 and all values 
of T and L y considered. We expect that w ~ L*, which 
we have checked explicitly by calculating w as the second 
moment of dm(y)/dy, showing that w decreases rapidly 
with Fq and eventually saturates at a finite value, see the 
inset in Fig. [2^i. Model II displays similar squeezing of 
the interfacial region, see Fig. [2}d. For small values of uo 
the effect is very similar to the boundary-driven case. In- 
creasing uo leads to very pronounced squeezing, such that 
the central section of the profile for uo = 1 resembles the 
mean-field solution in equilibrium! Even this case can be 
rescaled to the equilibrium profile, see Fig. [2}d. We find 
that the variation of w with transversal system size and 
field strength can be condensed into dependence on the 
scaling variable l/(L y uo s ) for any value of T considered; 
for T/T c = 0.75 data collapse is achieved for s = 0.3, 
see the inset of Fig. 2b. Similar quality of scaling can be 
obtained by including (T — T c )/T c in the scaling variable 
[7]; these results (not shown) indicate that in cases of 
weak drive a± = L*/L y increases (and hence the effec- 
tive confinement becomes weaker) with increasing T. 

We next investigate the interfacial structure on the 
two-body level and consider the spin-spin pair correla- 
tion function in steady-state G(x,y,y') = ((JiCTj}, where 
i = (0,2/) and j = {x,y'), focusing on the behaviour 
along the center, G\\(x) = G(x,L y /2,L y /2)] this should 
reveal most clearly interface- mediated correlations. The 
simulation results shown in Fig. |3^i indicate that shear 
induces a faster decay of G\\{x) with distance x than at 
equilibrium. For model I and for weak drive in model II, 
the behaviour at short and intermediate distances can be 
accurately described by rescaling the lateral coordinate, 
G\\(a\\x) w Gjj q (x). From fitting the data we find the ra- 
tio of the lateral (interfacial) corrrelation lengths in and 
out of equilibrium, ay = £||/£y q < 1, see Fig. |3| Unlike 
the behaviour on the one-body level, G\\(x) under strong 
drive in model II does not resemble the equilibrium mean- 
field result (obtained by calculating the inhomogeneous 
magnetization around a fixed spin in the interface), which 
displays significantly longer-ranged decay. 

We studied coarse-grained interfacial properties via the 
local position (height) of the interface h(x), obtained ei- 
ther following [17 or based on the (scaled) column mag- 
netization: h(x) = (2mb) _1 ^2 y <r(x, y)] both methods 
give consistent results. Fig. 3b displays results for the 
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Comparing our findings for C(x) to experimental re- 
sults for the height- height correlation function of a phase- 
separated colloid-polymer mixture, Fig. 3 of Ref. [5], re- 
veals a striking similarity in the reduction of the am- 
plitude of the correlation function under shear. This is 
consistent with the observation of a reduction of the in- 
terfacial roughness in both studies. However, Derks et 
al find an increase of the lateral correlation length under 
shear (Fig. 4b of [5]), as obtained from fitting the correla- 
tion function of the equilibrium capillary wave model to 
the steady state data. This finding is in contrast to the 
behaviour of our models. Despite the differences in the 
way the interface is localized (whether experimentally by 
gravity or here by walls) , the fundamental question arises 
as to which features of interfaces under shear are univer- 
sal. Clearly, dimensionality is expected to play a major 
role. The 3d Ising model has a rough phase above the 
(finite) roughening transition temperature. However, the 
roughness of the interface is established by a very differ- 
ent mechanism than in 2d, in that the interface throws 
out spikes [19], which are expected to hinder fluid flow 
to a much lesser degree than interface meandering in 2d. 
This should have important repercussions on the interfa- 
cial structure under drive. 

We thank D. Aarts, D. Derks, J. Eggers, R. Evans, 
A. Gambassi and R. Zia for useful discussions and the 
EPSRC for support. 



FIG. 3: (color online) Spin-spin pair correlation function 
G\\(x) at the center of the strip (a) and height-height cor- 
relation function C(x) (b) as a function of distance x for 
T/T c = 0.75, L x = 200, L y = 20. Shown are results in equi- 
librium, from mean-field theory (only for G\\(x) in a), from 
Kawasaki simulations with strong boundary-drive in model I 
(Fo = 50), the same rescaled as G\\(x/1A); Kawasaki simula- 
tion results for model II are shown for co = 0.025, 0.1 (only 
for C(x) in b) and 0.5, as well as the former case rescaled as 
G„(x/1.27). 

height-height correlation function C(x) = (h(x)h(O))', 
this is directly related to the spin-spin correlation func- 
tion: C(x) = (4m^) _1 J2y, y > G(x,y,y'). The intercept, 
C(0) = (h(0) 2 ), provides an alternative measure of the 
(squared) interfacial width; we find that (7(0) decreases 
under shear, very similarly to the behaviour of w 2 , ex- 
hibiting a very strong effect for large values of uj. C{x) 
decays faster with x than in equilibrium, in agreement 
with the behaviour of the spin-spin correlation function. 
Hence from the scaling of m(y) and G\\(x), we expect 
a~]_ 2 C(a || x) « Ceq(x), in correspondence with Weeks scal- 
ing in equilibrium [18]: C eq (x) ~ w 2 C(x/^ q ), where C 
is the scaling function, and w ~ L y . Indeed both for 
strong boundary drive and weak bulk drive data collapse 
is achieved using the same values for ay and a± as ob- 
tained above. 
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